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MUTUAL INDUCTANCE OF ANY TWO CIRCLES 

By Chester Snow 



ABSTRACT 

A formal expansion is derived for the magnetic flux through any circle. A 
similar expansion is obtained for the magnetic potential due to a unit current 
in another circle. Combining the two gives a formal expansion for the mutual 
inductance of the two circles. Application is made to two cases (1) where the 
circles are parallel, (2) where their axes intersect. 



CONTENTS 

Page 

I. Introduction 531 

II. Symbolic expression for magnetic flux through a circle 532 

III. Symbolic expression for the magnetic potential of a circular current- _ 535 

IV. Combination of II and III for the mutual inductance of the two 

circles 536 

1. Case of parallel circles 536 

2. Special case where the axes of the two circles intersect 538 

(a) The axes intersect at center of one circle 540 

(6) Circles coaxial 541 

(c) Circles concentric, but not necessarily coaxial 541 

V. Convergence of the expansion 541 

I. INTRODUCTION 

In certain absolute measurements it is necessary to compute with 
precision the electromagnetic forces exerted between two current- 
carrying solenoids, or the current induced in one by changes of that 
in the other. Both problems require a formula for the mutual 
inductance of the two, which is generally to be derived by integrating 
the expression for the mutual inductance between two circles. The 
fundamental importance of the two circular elements is evident. 
When the circles are not only parallel but coaxial, their mutual 
inductance is given exactly by Maxwell's formula. 

In 1916 Butterworth ! derived a formula covering (partially) the 
cas.e of two parallel circles. In the present paper is presented a general 
formal expression for any two circles which leads readily to Butter- 
worth's formula as a special case, and also to a fairly simple expansion 
in zonal harmonics for the mutual inductance of two circles whose 

i S. Butterworth, Phil. Mag., p. 448; May 31, 1916. 
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axes intersect at a finite point (equation (23) below). From it may be 
derived simple formulas for the mutual inductance and torque between 
two solenoids whose axes intersect, one lying entirely within the other, 
as in certain types of current balance. 

The starting point of these developments is a formal expression 
here obtained for the magnetic flux through a circle of radius a 2 , the 
field being the negative gradient of a scalar potential 12 (x 2} y 2 , z 2 ) at 
the point P 2 (x 2 , y 2 , z 2 ), which is the center of the circle. It is 

M= - f f D n2 QdS = -2ra 2 J x (a 2 D Q2 ) 12 (x 2 , y 2 , z 2 ) 

where J\ is BessePs function of the first order, the symbol D n2 des- 
ignating the space derivative in the direction normal to the plane 
of the circle — at the center. 
Thus 

oo / a \ 2k+l 

Ji (a 2 D Q2 ) o - Y] (-Wyi) D D 

£t Hf+1)! 
where 

D n2 Q = (l 2 D X2 + m 2 D y2 + n 2 D Z2 ) Q 

the direction cosines of the normal to the circle being l 2 , ra 2 , n 2 . 

If the field 12 (x 2) y 2 , z 2 ) is that produced by a unit current in a circle 
of radius a x center at Pi {x% } y\, Zi), it is then shown that 

12 (x 2 , y 2y z 2 ) - 2*0* Ji (0|Z? B1 ) - 

where r is the distance between centers. The mutual inductance 
between the two circles is then given (formally) by 

M= -^axOiJi (a 2 D n2 ) J x (aiZ? ni ) - 
By obtaining expansions for - which are suitable from the point of 

view of performing the infinite series of differentiations here indi- 
cated, a variety of series formulas may be obtained. 

II. SYMBOLIC EXPRESSION FOR THE MAGNETIC FLUX 
THROUGH A CIRCLE 

Let 12 (as, y, z) be a scalar point function which may be expanded by 
Taylor's (three-dimentional) theorem about the point x , y , z . This 
expansion may be written symbolically if we let 

x = x +x', y=y +y', z = z + z' 

12 (x, y, z) = 12 (x + x', y Q + y',z + z') - e ^^ Y '^ +z ^ 12 (x , y oy z ) (1) 
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Consider a circle of radius a whose center is at x , y oy 2 and whose 
plane is perpendicular to the z axis. Let u denote the surface inte- 
gral of £2 over this circle. Then 



u= | dx' I dy f (x + x', y + y'y z + z') 

-V* 



J a /»Va J -x ,J 

fo' I %' Q ( 
a «/_ A / a 2_ x 'j 

/*a /»Va 2 "^V2 

I cfc' I (ft/' e x ' D "o+ y 'D yo ft (a- y 0j 2o ) 

•/-a t/Va 2 -x /J 



(2) 



"V8 

In this expression the symbolic operator e*' Dl ° +y ' Dj, ° Q (x , y , z ) 
is merely the abbreviated form of the series of operations indicated 

by 

e *<D, 0+ y'D yo Q (a . o; yo) Zo) = ^ji?^o+^AJ_ n {Xo> ^ z<>) (3) 

n=o 

where each term (x'D^ + y'DyJ 111 to (z , 2/ , 2o) is to be expanded for- 
mally by the binominal theorem. The result of using this formal 
expansion equation (3) in the integral, equation (2), may be found by 
noting that the symbols D Xo and D Yo behave like constants as far as 
integration with respect to x f and y' is concerned. It may be written 
out from a consideration of the formal analogy between it and the 
integral u where 

J a /"Va 5 " 11 ?! 

dx f I dy' e x ' A +y'B a o (4) 

a J _ ^U^ 

where A, B, and & do not involve x' or y' and are, therefore, con- 
stants as far as this integration is concerned. This integral may be 
readily evaluated by choosing a new pair of rectangular axes x" , y" 
in the same plane as before and with the same center such that 

n , a , ; B 



-yjA 2 +B 2 y -y/A 2 + B 2 
In terms of these axes the expression for u Q is 



dx" I dy" e*' f JA*+B> = 2a 2 tt I s 

a t/_7 a j- x 'i t/ o 

OO (*TC 

= 2ou£] fe?SB35E I cos" sin 2 dd9 

n = o t/ o 
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Now I cos n 9 sin 2 ddd = if n is odd 



■: 



— if n is even; n = 2s; s = 0, 1, 2, 3- 



2 2s + 1 s!(s + l)! 
Hence 

2 >,V 2 /_ _ J 1 (m^ + m a ... 

= I I g x ' A +y' B & dS integrated over the circular area. J\ 

is Bessel's function of order one. 

A comparison of this equation with the symbolic equation (2) en- 
nables the latter to be put in the form 



f f fi(0Vs„+^„) s 






) 

(6) 



The meaning of ^D 2 Xo + D 2 yo is merely an operator whose repetition 
is equivalent to the operator D\ + D 2 yo , and since none but positive 
even "powers" of this operator occur here it need not be further 
defined. 

The expression (6) assumes a particular interest in case the func- 
tion ft (x, y, z) satisfies Laplace's equation for in that case 

(D\ o + D\ o yn(x , Vo , s ) -(- l)*D\Mxo, 2/c to) (7) 

and (6) becomes 

I [ MS - wa'f) ( " 1)8 ( V 2 ) Qfo, y , z ) (8) 

J J s=o «!(k+i)l 

The surface integral of any harmonic point function £2, over any 
circle of radius a, is thus expressed in terms of the values which 
are assumed at its center by the function and its normal derivatives 
of even order. If 12 (x, y, z) is harmonic, then so is D Zo Q. Hence, 
if we replace £2 by D Xo 2 it becomes 

z>, o o^=2^V(-i) a V^-; aOc,y.,go) 

J J fc? «!(«+l)! (9) 

= 2TaJi(aD Zo )tt(x , y ol s ) 
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If fi represents the magnetic, scalar, potential whose negative gra- 
dient is a magnetic field, then the magnetic flux through the circle 
may be written in the symbolic form 

M=-ffD n ndS= - 27ra 2 J 1 (a 2 D Q2 Mx 2 , y 2 , z 2 ) (10) 

where a 2 is the radius of the circle, P 2 (x 2 , y 2 , z 2 ) its center and D n2 Q 
represents the value at P 2 of the space derivative of 12 in the direc- 
tion of the positive normal to the plane of the circle, which is such 
that the positive direction of the current encircles it right-handedly. 
It is evident that M is harmonic in the variables x 2 , y 2 z 2 , if the orien- 
tation of the circle is held constant; that is, (D 2 x2 + D 2 y2 -\- D 2 X2 )M >-=Q. 

III. SYMBOLIC EXPRESSION FOR THE MAGNETIC POTEN- 
TIAL OF A CIRCULAR CURRENT 

If there is a unit current in another circle of radius a u whose plane 
is in the yz plane and whose positive normal is the positive x axis so 
that its center is at the origin of coordinates, then the (scalar) magnetic 
potential Q, (x, y, z) due to it at any point P(x f y 1 z) may be found 
from the fact that it is harmonic and is a function of the two cylin- 
drical coordinates only, namely, x and p^-\/y 2 + z 2 } and reduces 
when x = + to if p^>ai and to 2 w if p<C^i. (This is evident since Q 
is the solid angle subtended at P by the positive area of the circle.) 
From the known discontinuous function defined by the integral 

CO 

2Tra l fJ i (a l s)J o (ps)ds = if p>a x 
= 2x if p<^ai 

and from the fact that a function of x and p only, which is a particu- 
lar solution of Laplace's equation, is e^ BX J (p s), it is readily seen that 



Q(x, p) =2wa 1% fe BX J 1 (ais)J Q (ps)ds if z>0 



CO 

= —2iraiJ t e ax Ji(a 1 s)J (ps)ds if #<0 



(12) 



where 

P=V2/ 2 + * 2 (13) 

The expression for ft may be given another form by using the formula 

PV iW Jo (s V^ 2 ) ds= \ (14) 

Jo -yjx 2 + y 2 + z 2 

and by noting that, since 

{-a x D x ) n e-**={a x 8y e~* x 
-J x (oiDJ e~ Si = J l (-a x D x ) *-»-*-«/, (a^) 
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and consequently 

e~ sx J x (a x s) J (ps) ds=-J x (a x B x ) e~ sx J (ps) ds 

Jo Jo 

•y/x 2 +y 2 + Z 2 

Consequently, the two expressions (12) for ft may be put in either 
of the following forms, each of which holds for both positive and 
negative values of x, namely 

fl (x, p) = -2ira l J x (a x D x ) , j— 

^x 2 -\ry 2 + z 2 

ro, (15) 

= -2ira x $\ (a x D x ) J «- 8 W J (ps) ds 

Remembering the special relation of this circle to the x axis, it is easy 
to write out from a consideration of equation (15) the expression for 
the value of the magnetic potential 12 at a point P 2 (x 2 , y 2 , 2 2 ) due to unit 
current in a circle of radius a x with center at P x (x X) y x , z x ), whose plane 
has the positive normal n x . It is 

(*,, y 2 , z 2 ) = 2^ J> ffgjj Vte-xO'+dh-yJ'+fa-eO' ° 6) 

where a x D Ul represents the directional derivative in the direction n X) 
with respect to the variables x x , y Xj z x \ that is, if l h m x , n x are the direction 
cosines of the positive normal to this circle 

a x D ni = a x (l x D Xl + m x D yi + n x D zl ) 

IV. COMBINATION OF II AND III FOR THE MUTUAL IN- 
DUCTANCE OF TWO CIRCLES 

Equation (16) is purley formal, like equation (10), to which it 
bears a certain resemblance. Combining the two, leads to the 
following symbolic formula for the mutual inductance of any two 
circles of radii a x and a 2 with positive normals n x and n 2 . 

M= -4T 2 a x a 2 J x (a x D Ql ) J x (a 2 D n2 ) , - (17) 

■>j(x 2 -x x ) 2 +(y 2 -y x ) 2 +(z 2 -z x ) 2 

where P x (x Xy y X) z x ) and P 2 (x 2) y 2 , z 2 ) are the centers of the circles, respec- 
tively, and n Xy n 2 are unit vectors normal to their positive sides. 

1. CASE OF PARALLEL CIRCLES 

In the special case where the circles are parallel, the first one may 
be taken in the yz plane with its center at the origin, the center of the 
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second in the xy plane. Denoting the coordinates of the second 
circle by xy (or polar coordinate r } 6), the formula (17) gives 

M=4:7r 2 a l a 2 Ji (a x D x ) J x (a 2 D x ) - 

/*oo 

= 4tt 2 a x a 2 I e~ slx| J Q (ys) J x (a x s) J x (a 2 s) ds 
« F(-n-n+l,2,p) 

n=l 

where F denotes the hypergeometric function. But 



1 (27i)!„ . M _2 2n " 



!r ( n ^) 



so that 



D 2 \ f rjg? P*» (cos*) =-j= - ^ 7 P ia (cos©) 

oo 
v n = l 

(n— 1)! \ ' ' ' cr,/ 



(18) 



which converges for large values of r. This is identical with a formula 
obtained by Butterworth. 2 

Formula (18) is valid when r>a l + a 2 . 

In the other extreme, where r<a { — a 2 , one finds 



r _^yj, r :Jb§ 



M- 

(18') 
^ + i. + |2,g-;)(0P 2n (cos.) 

This may be derived from (17) by writing the latter in the form 

M= —Air 2 aia 2 J (yD x ) J x (a 2 D x ) J e~ BX J x (a x s) ds 

= 47T 2 a 2 Jo (yDx) J i (a 2 D x ) , 2 

-yd i X 

by using the binomial expansion 






2 See footnote 1, p. 531. 
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together with the operational expansion 

J (yDMa 2 DJ = -fZj (- !)■- , ^/i) , 

2. SPECIAL CASE WHERE THE AXES OF THE TWO CIRCLES 

INTERSECT 

A considerably more general case (but not the most general) is 
that in which the axes of the two circles intersect. Taking the plane 
of these two axes as the xy plane, with the origin at the point of inter- 
section of the two axes, the center of the first circle being on the x 
axis with abscissa Xi (xi may be taken as never negative without loss 
of generality), and that of the second circle having the plane polar 
coordinates r, 0, then (17) becomes (see fig. 1) 

M = - 2xa 2 Ji (a 2 D T ) Q (r, 0) 

- - 27^! (a 2 D T ) 2ira l J l (a^) , - ---—-. 

V a i — % x i r cos + r 
There is no loss of generality in taking for circle No. 1 that one of the 
two whose center is most distant from the origin. Hence, we may 
use the expansion 

oo 

1 \"1 r 8 

/ 2 o /> i Iz = / i~ i+i Pb ( cos 0) since £i>r 

V» i — 2x 2 r cos + r 2 iLJxi B+1 K • ^ 



s=o 



where P B is Legendre's coefficient. Applying to this expansion the 
symbolic formula (16) one finds 

v.,(-««(ir 

_o_^, / J \zLL 71 2k+l _ 

OO 00 

_ _t^ VI VI (-l)*r(2fr + s + 2) / aj_Wr V 

x ! ,ZjZjr(i+i)r(it+2)r(s+i) ^2*1/ W ' l ; 

k=o s=o 

-pS(j,)' p -<-«»V^TT)]- 
Lfe 1 r(*+i)r(fc+2) v Xl/ J 

xa 2 iV^/r\ s , , 1Nto , * »/*+2 s + 3 „ a 2 ,\ 
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00 

sin 2 ai][](£) S ( s+ W ^ (cos 0) f(^p ^=p 2, sin 2 a,) 



= — 7r sm' 

s=o 



where 



tan «i =— * and r 2 i = x 2 x 4- a 2 i 

#1 



Now by Gauss's transformation one finds that 

jp/2 + 8 1 — 8 . 2 \ 2(l-C0S«i) •/. n 

F (-2-' X' 2>Smai r *n'«i lf5 = 

_ 2P^B« 1 ) if (21) 

s(s-fl) ' ' ' v ' 

where P' 8 denotes the derivative of P 9 with respect to its argument. 
Using (21) in the preceding expression for Q gives 



Q(r, 6) = - 2tt(1 - cosaO - 27rsin 2 «iJJ (r")' 



r V P 8 (cos 0)P' 8 (cos aO (22) 



a formula which holds when the origin is any point on the axis, 
to the left of the circle if r<^r x . Substituting this value of £2 in the 
symbolic formula (19) gives 



v 2k+l 2k+l 



/ n \2k+l 21 

u a 2 • 2 VI V V V V a P B (cos fl)P B (cos ad 
M=^a 2 sm*a 1 2j £!(fc+l)! 2j ~^7 



k=o s=l 

v 2k+l 



: 47T 2 a 2 sin 2 a 1 

k=o s=l 



or 



- » (-i)^|) 2k+I (s -i)! 
2j 2j *!(*+l)!r(«-2*) (rj P » (cos *) P '-( cos «0 

k=o s=l 

OO CD f— 1U/^?Y 

27rWsin 2 a i Yl^ry D ,yi *■ 1; ^ 

s=l k=o 

oo 

M ^flVj yi/r 2 V P B (cos flP'.fcos « 1 )P / .(cos a 2 ) , . 

7yr 2 ZJvn/ 5(5 + 1) l ' 



where 



r 2 ! = x 2 i+a 2 i, tan «*«— l 



a, 

(See fig. 1.) 



r 2 2 = r 2 + a 2 2 , tan a 2 = — 



«i and a 2 are positive acute angles. 
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The formula (23) gives the mutual inductance between any two 
circles of radii a t and a 2 whose axes intersect. The circle No. 1 is 
chosen as that one of the two whose center is most distant from the 
origin. The origin is the point of intersection of the two axes and 
the angle between these axes and may be anywhere in the range 
from zero to it. The origin is taken as lying to the left of the circle 
No. 1 whose center is on the x axis at x = x x . The distance of the 
center of the second circle from the origin is r. (See fig. 1 where 
r = 00 2 .) 

(a) The Axes Intersect at Center of One Circle. — A special 
case of (23) is that in which the center of the second circle 2 coincides 




Fig. 1. — A principal section of two circles whose axes intersect 



IT 

with (fig. 1) so that r 2 becomes a 2 and a 2 becomes =»• In this case 
P' s (cos a 2 ) vanishes when s is an even integer, and if s is an odd 



integer, say 2n+ 1, reduces to 
so that (23) reduces to 



VA / Til 



<*4) 



M= 2tt! sin «i V^2j( "" 1 ) \7) (n+l)i P2a +^ cos ^ pf ^ +lCOSa ^ 

0=0 

(23a) 
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(b) Circles Coaxial. — Specializing this still further by making 
the circles coaxial (so that £i=cos 6*=1) gives 



M=2*th sin ^V / 7ry^(-l) n ^V n+2 r ( 7l + 27 
td W 

when 



(n+1)! 



PWcos a x ) (23b) 



a 2 



<1 



This is easily shown to be reducible to Maxwell's formula involv- 
ing elliptic integrals. 

(c) Circles Concentric, but not Necessarily Coaxial. — 
Specializing (23a) in another direction, by letting the circles become 
concentric (a 2 being less than a x ) y gives 



M: 



■ ^Wa)- ^ + gX^ + |) p M(cOB d) (23c) 



Thus, by considerations of continuity, we have arrived by trans- 
formations at a formula which is valid, where the original assump- 
tion with which we started is not satisfied, namely, that the field of 
a 2 is expressible by Taylor's theorem over a sphere of radius a x . 

V. CONVERGENCE OF THE EXPANSION 

Returning to a consideration of the general formula (23), it is 
readily seen that this series converges when r 2 <Cr u so that by inter- 
changing r 2 and r x a formula is obtained for the case when r 2 ^>r x . 
The proof of the absolute convergence of this series may be based 
upon the fact that P 8 is never greater than unity when its argument 
is the cosine of a real angle and upon the fact that 



P' 8 (cosa): 
which shows that 



g(g+l) 

•4 



P e _ (cos a) — P B+ i(cos a) 
2 sin 2 a 



P B (cos fl)P' s (cos cQP'^cos a 2 ) 
8(8+1) 




r Ps-^cos a x ) - P 8+ i(cos aj l pVxfcos a 2 ) - P a +i(cos a 2 ) l 
L 2sin 2 a! JL 2sin 2 o: 2 J 



r 2 ! i 



: sin 2 <*! sin 2 a 2 a 2 id 2 2 
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Hence, if r 2 <jr l 

CO 

*-J\rJ r x -r 2 

s = l 

♦ 

No discussion of the convergence range' of the symbolic formula 
is possible without first choosing the type of coordinates suitable for 
performing the series of differentiations. The convergence range 
usually makes itself evident in the process of performing the 
expansions. 

Washington, May 19, 1928. 



